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Abstract: 

In this work, we present an alternative subtraction scheme for next-to-leading order 
QCD calculations. The subtraction terms are derived from the splitting functions of an 
improved parton shower with quantum interferences, and our scheme exhibits a reduced 
number of momentum mappings in the real emission subtraction terms and facilitates 
the mapping of NLO calculations with the shower. While simple cases with maximally 
Q-i two partons in the final state have already been presented in a previous publication, 

we here complete the discussion of the massless scheme by providing the final state real 
emission and integrated subtraction terms for processes with any number of partons 
in the final state. We apply our scheme to three jet production at lepton colliders at 
next-to-leading order and present results for the differential C-parameter distribution. 



1 Introduction 



With the start of data taking at the LHC in 2009 and its more than successful physics 
program since then, particle physics has entered an exciting era. Major tasks of the 
LHC experiments are the accurate measurement of the parameters of the Standard 
Model (SM) of particle physics, as well as the search for physics beyond the SM (BSM). 
For both, a precise understanding of the SM signals and background processes in an 
hadronic environment are crucial. These processes are mainly governed by strong in- 
teractions, where leading order (LO) calculations can exhibit uncertainties up to 100 % 
(c.f. pQ for a recent review); therefore, for a correct theoretical prediction of these 
processes at least next-to-leading order (NLO) corrections need to be taken into ac- 
count. Furthermore, it is generally not sufficient to apply an overall NLO rescaling 
K-factor to the leading order predictions, as NLO corrections can vary widely for dif- 
ferent regions of phase space. More accurate predictions therefore call for the inclusion 
of these NLO calculations in Monte Carlo Event generators, which provide predic- 
tions for fully differential corrections to the LO process. Many such generators exist 
at parton level [2HT8] , and in recent years a lot of progress has equally been made to 
automatize the matching of these processes with parton showers in the Powheg [THH3S] 
and (a)MC@NLO [22JII7HS5] framework^- 

In this paper, we present the generic extension of an improved subtraction scheme 
[58H60] . which facilitates the inclusion of infrared (IR) NLO divergences originating 
from different phase space contributions in Monte Carlo Event generators. These di- 
vergences arise whenever internal loop momenta approach zero or particles become 
collinear, and are known to cancel in any fixed order in perturbation theory [HUE2]- 
However, these cancellations occur in the combined sum of virtual and real emission 
contributions, and therefore originate from phase spaces with a different number of 
particles in the final state. In analytic or semi-numerical calculations, the singulari- 
ties can be parametrized by an infinitesimal regulator; in the sum of real and virtual 
contributions, these regulators can then be set to zero to obtain a completely finite 
prediction. In numerical implementations, however, the inclusion of infinitesimal reg- 
ulators can easily lead to numerical instabilities. Subtraction methods [63HT3] circum- 
vent this problem by introducing local counterterms which mimic the behaviour of the 
real emission matrix elements in the singular limits. The integrated counterparts of 
these terms are then added to the virtual contributions, where again an infinitesimal 
regulator is used to parametrize the singularities. Then, the higher order contributions 
in both phase space integrations are finite, and the regulator can be set to zero. In 
recent years, many of these schemes have been made available on a (semi) automated 
level $Nn\UM3\. 

While the behaviour of the subtraction terms in the singular limit is determined by 
factorization [75H50] . the finite parts of the local counterterms as well as the mapping 
prescription between real emission and leading order phase space kinematics in the 
subtraction terms can differ. Unfortunately, standard schemes [65l[70] suffer from a 
rapidly rising number of momentum mappings, which scales like N 3 for a leading order 

1 Recent reviews on this can be found in [551157] , 



1 



2 — > N process. Therefore, increasing the number of final state particles leads to a 
rapidly rising number of reevaluations of the Born matrix element. In [q*%HBT)] . we there- 
fore proposed a new subtraction scheme with a modified momentum mapping [STW85] . 
where the number of momentum mappings scales as ~ iV 2 . The momentum mappings 
are constructed such that they take the whole remaining event as a spectator, and 
the subtraction terms are derived from the splitting functions in an improved parton 
shower [5TH55] . In [58HBT)] . we presented the scheme for the simplest cases with max- 
imally two partons in the final state^. In the present work, we extend the scheme to 
cases with an arbitrary number of massless particles in the final state. We also provide 
the helicity dependent squared splitting functions for splittings where the mother par- 
ton is a gluon. We validate our scheme by applying it to three-jet production at NLO 
at lepton colliders, obtaining complete agreement with the Catani Seymour scheme. 

This paper is organized as follows. In Section 2, we briefly review the generic setup for 
subtraction schemes. In Section 3, we review the ingredients of the new scheme and 
present the generalized results for the integrated subtraction terms for an arbitrary 
number of final state massless partons. We discuss the application of our scheme 
to three-jet production at lepton colliders in Section 4. Conclusion and outlook are 
presented in Section 5. The appendix contains a summary of the final state splitting 
functions [SI] used as subtraction terms, a generic parametrization of four-parton phase 
space, and the collinear subtraction terms for processes with incoming hadrons. 

2 General structure of NLO cross sections and sub- 
traction schemes 

In this section, we briefly review the general subtraction procedure for calculating NLO 
cross sections at lepton and hadron colliders. We start with a generic cross section at 



where a should be specified by the respective jet function as discussed below, and 
da B , da R , and da v are the Born, real emission, and virtual contribution respectively. 
We here consider processes with m particles in the Born-contribution and m + 1 partons 
in the real emission terms. After UV-renormalisation, the virtual and real-emission 
cross sections each contain infrared and collinear singularities. These cancel in the 
sum of virtual and real contributions [6TJ[62], but the individual pieces are divergent 
and can therefore not be integrated numerically in four dimensions. 
Subtraction schemes consist of local counterterms which match the behaviour of the 
real-emission matrix element in the soft and collinear regions, and their integrated 
counterparts. Subtracting these counterterms from the real-emission matrix elements 
and adding back the integrated counterparts to the virtual contribution results in 

2 Some results for the generic scheme were already presented in [58| . 
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Figure 1: Soft /collinear factorisation: when the partons i and j become collinear 
and/or parton j becomes soft, the (m + l)-parton matrix element factorises into a sum 
over m-parton matrix elements times a singular factor vg. 



finite integrands for both the virtual correction (m-particle phase space) and the real 
contribution ((m+ l)-particle phase space): 

da NLO = [da R - da A ] + [da A + da v ] . (2) 

The construction of the local counterterms, collectively denoted by da A , relies on the 
factorisation of the real-emission matrix element in the singular (i.e. soft and collinear) 
limits (Fig. [T]) [75H5U] . and we symbolically write 



\M m+1 (p)\ 2 — > V,®\M m {p)\\ (3) 

where T>i are the dipoles containing the respective singularity structure, and the sym- 
bol ® denotes a correct convolution in color, spin, and flavour space, ft/ p represent 
momenta in [m + 1)/ m-parton phase space, respectively. As |A^ m +i| 2 and |.M m | 2 live 
in different phase spaces, a mapping of their momenta needs to be introduced, which 
is defined by a mapping function F map according to 

P = F mSLP (p). (4) 



T>i and its one-parton integrated counterpart Ve are related by 

Vt= f d£ p V t , (5) 



where d£ p is an unresolved one parton integration measure. 



In summary, any subtraction scheme needs to fulfill the following requirements: 

• The dipole subtraction terms T>i must match the behaviour of the real emission 
matrix element in each soft and collinear region, and lead to correct IR poles 
when carrying out the analytical integration over the one parton phase space in a 
suitable regularization scheme that are necessary to cancel the soft singularities 
in the virtual (one-loop) matrix element, 
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• The mapping function F map guarantees total energy momentum conservation 
as well as the on-shell condition for all external particles before and after the 
mapping. 

Integrating Eqn. (J2J) over phase space using dimensional regularization [84,85], where 
D = 4 — 2 e, then yields 



a 



NLO 



f [da R - da A ] + I da A + I 

Jm+l Jm+l J rr 



da 



V 



finite 



finite 



m+1 



[daf =Q - da A =0 ] 



da A + da 1 



(6) 



J £=0 



Both integrands are now finite: the integration in (m + 1) particle phase space can 
safely be performed in D = 4 dimensions, as the singular regions are regularized by 
the respective counterterms. In the m parton phase space, the sum of the integrated 
dipole contribution and the virtual correction does not contain any further poles, so 
that we can set e — 0. Then, all integrations can be performed numerically The 
explicit expressions of the cross section a for m and [m + 1) particle phase space 
contributions at NLO are 



da a + da v + 







f dPS m 









\M m \ 2 + \M 



m I one-loop 



[ [da R - da A ] = [ dPS m+l 

Jm+l J 



\M 



m+1 1 



- D t (8) \M, 



(7) 



where in this symbolic notation j dPS includes all flux and symmetry factors, and with 
|-M m | 2 , (A^m+il 2 and |A^ m | 2 ne _ loop being the squared LO matrix element, the squared 
real emission matrix element and the interference term, respectively. In Eqn. (JTj), the 
sum runs over all local counterterms needed to match the complete singularity structure 
of the real emission contribution, and convolution with jet functions then ensures the 
collinear and infrared safety of the Born-level contribution. The insertion operator 1(e) 
is then defined on a cross section level according to 




da J 



da 1 



Ke) 



where the symbol <%> again denotes a proper convolution in spin, color, and phase space. 
The generalization of this for processes with initial state hadrons has already been pre- 
sented in [59]; for completeness, we repeat the argument in Appendix IBl 
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Observable-dependent formulation of the subtraction method 

The jet observables should be well defined such that the leading order cross sections 
are infrared and collinear safe. The jet cross sections are defined as 



a 



LO 



a 



NLO 

J 



I dPS m (pi, ■ ■ ■ ,p m ) \M m {pi, ■ ■ ■ j Pm ) | (Pl> • • * >Pm), 

dPS m+ i(pi, ■ ■ ■ ,p m +l) \M m +l(pi, ■ ■ ■ ,Pm+l)| 2 FjT (Pl> " " " >Pm+l) 



+ J dPS m (px,--- ,p m ) \M m (p 



I,-' ,P 



F j m) (Pir-- ,Pm)- 



m ) lone-loop 



In general, the jet function may contain 9 functions (which define cuts and correspond- 
ing cross sections) and 5 functions (which define differential cross sections). For an 
infrared finite jet function, we require that 

f j (Pi,--- ,Pj = A?,-" ,Pm+i) -»■ Fj m) {Pi, ■ ■ ■ ,Pm+i), if A^O 

(PU~,Pi, ■;Pj,-;Pm+l) ~+ F^ (Pi , ■■ , P, ■ ; Pm+l) 

if Pi ->■ zp, pj ->■ (1 - z)p, 

F$ m) (l>ir-- ,Pm) ->• 0, if Pi-p^O. (9) 

The last condition of Eqn. (j^J) corresponds to an infrared safe definition of the Born- 
level observable, while the first two conditions guarantee infrared and collinear safety 
of the observables and can be summarized to 



in the singular limits. 



(10) 



We then have 



da B + da v -l / da J 



dPS n 



\M m \ 2 + \M m \ 2 onc _ loop + ^V,® \M 7 



f!T\p), 



m+1 



[da R - da A ] 



dPSm+i 



\M m+1 \ 2 if 1 +1 \p) -£>® IM^F^ip) 



'12) 



(13) 



where in the integrated subtraction term the momenta p are derived from p using the 
respective momentum mapping. 



3 Alternative subtraction scheme: setup 

In this section, we will first review the setup of our scheme as well as the mapping and 
respective integration measures which have already been presented in [59l[8T]. In our 
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scheme, the NLO subtraction terms are derived from the splitting functions introduced 
in a parton shower context [8"TH8"3] , and the m+ltom phase space mappings needed 
correspond to the inverse of the respective shower m to m + 1 mappings. In the 
following, we will denote the m+ 1 phase space four- vectors by pi,pz, ■■■ and m phase 
space four-vectors by pi, P2, .... Inm + 1 phase space, the four-momenta of the emitter, 
emitted particle, and spectator are denoted p£, pj, and f>k respectively. Note that here 
the spectator only needs to be specified if pj denotes a gluon, as we use the whole 
remaining event as a spectator in the sense of momentum redistribution for both initial 
and final state mappings. We here restrict our expressions to subtractions on the parton 
level and to massless partons. 



3.1 Splitting functions 

We start with a description of the matrix element factorization in the soft and collinear 
limits, following the notation in [81] , where the QCD scattering amplitude for m + 1 
partons is a given as a vector in (colour (g> spin) space, 

\M({pJ} m+ l)). (14) 

In the singular limits, the amplitude | Aie({p, f} m +i)) can be factorized into a splitting 
operator times the m-parton matrix element 

| M t ({p, /WO) = Aih ^fe + fj) V}{{p, /WO I M({p, /} m )>, (15) 

where the index I labels the emitter/ mother parton in the (m + 1)/ m particle phase 
space. V^({p, /} TO +0 i s an operator acting on the spin part of the (colour g) spin) 
space, while t\(f t — > fi + fj) is an operator acting on the colour part of the (colour £g> 
spin) space. The Born amplitude for producing m partons is evaluated at momenta 
and flavours {p, f} m determined from {p, f} m +i according to the respective momentum 
mappings. The spin-dependent splitting operator can be described in the spin space 

({sWl I V?({p,/Wl) I Wm>- (16) 

If we take Eqn. (fT6|) to be diagonal, we can define the splitting functions vg according 
to 

({s} m+ i\V e \{p,f} m+1 )\{s} m )= j 1 j 5§ n;Sn J v e ({p,f} m+1 ,Sj,Si,si). (17) 

Explicit forms for the splitting functions vg have been presented in [ST]. For the con- 
struction of the subtraction terms, we consider the approximation for the squared 
matrix element in the singular limits 

E (M e ({p, / WilMHiP, /WO ~ E v\ve (M{{pe, f} m )\M({pe, f} m )). 

For the direct splitting function, where I = £', we obtain 

W u = v j = vg({p, f} m+1 ,Sj,s e ,se)Vg({p, f} m+1 ,Sj,s e ,se), (18) 
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which, after summing over the daughter parton spins and averaging over the mother 
parton spins, leads to the spin averaged splitting functions Wu as subtraction terms. 
If the mother parton is a gluon, the Born-type matrix element might have an explicit 
dependence on the gluons polarization; in this case, we need to use 

(u\WuW) (19) 

in the real-emission subtraction terms, where u, v' are the polarisation indices of the 
m-parton phase space gluon. If the spin correlation tensor defined by Eqn. f lT§j) 
is perpendicular to pg, the angular correlations vanish after the integration over the 
unresolved particles phase space and the integral over Wu still provides the correct 
integrated counterterm |65j. For the collinear terms, the colour factors can easily be 
obtained [8T] : 

{C F (ft, ft) = (q,g),(g,q), 
c A (fijj) = (9,9), 
t r (fi,f S ) = (q,q). 

For soft gluon emissions, we also have to consider terms for which I 7^ £', which we 
will describe below. 



3.1.1 Eikonal factor 

When a gluon with four-vector pj becomes soft, or soft and collinear with p£, the 
splitting amplitude vi defined in Eqn. fTl7|) can be replaced by the eikonal approximation 
for pj — > 

vf\{P, /Wi, Sj, s e , at) = V4^Ts8 tM e fc'fi'Q)*-fr ( 20 ) 

Pj ' Pi 

where s(pj, Sj, Q) denotes the polarization vector of the emitted gluon with spin Sj. Q 
denotes the total momentum of the (m+1) phase space event and is used as a gauge 
vector. The eikonal approximation of the spin-averaged splitting functions Wu is then 

W eik Avr n Pf D (PjiQ) 'Pi / 9 ,X 

W u = 4ira s — , (21) 

VP 3 ■ Pi) 

where flavour-dependent averaging factors are already taken into account. The trans- 
verse projection tensor is given by 

Pj ■ Q {Pj ■ Q) 2 

It will be convenient to define a dimensionless function F: 

F = P ^W U . 
47T a s 

We then have 

7-1 _ Pe- Pj TTTcik pe ■ D(p j} Q)-p t 2p e ■ Q Q 2 p t ■ pj 

-Teik = - A yy ee ~ ' — ' — ~ " — 7\ T~~ — TTvT- 

Aixa s Pi -Pi Pj-Q {Pj-Qr 

The eikonal factor, in combination with the interference terms, is then used to construct 
dipole partitioning functions. 
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M t ({p,f) 



(M k {{p, f} m +i) 




I m + 1 



Figure 2: Soft diagram: parton j is emitted from parton I in the scattering amplitude 
and parton j is emitted from parton k in the complex-conjugate scattering amplitude. 



3.1.2 Soft splitting functions 

For soft gluon emission, we also need to take interference diagrams between different 
emitters into account. This means the emitted parton j can be emitted from emitter 
£ in the amplitude and parton j can also be emitted from a different emitter k in the 
complex-conjugate amplitude (Fig. [2]). The interference splitting function is then given 
by 

W tk ~ v e ({p, f} m+ i, Sj, s £ , Si) v k ({p, f} m+ i, Sj, s fc , s k )* 8 htSl 5g k , Sk . (22) 

The splitting function Eqn. (1221) contains a singularity when the emitted gluon j is soft; 
however when gluon j is collinear with parton £ or k, it does not contribute a leading 
singularity. In the special case that pj is soft, or possibly soft and collinear with p£, we 
can use: 

W a ~ v c e lk ({p, /} m+ i, fy, h, s^ v% k ({p, f} m+1 , sj, s k , s k )* 5 h)St 5s k , Sk . 

Note that this term only contributes if particle j is a gluon. In this prescription, there 
is an ambiguity in the allocation of the singularities, which can be distributed with the 
help of dipole partitioning functions; for completeness, we here repeat the argument 
in [81j[82]. The complete sum over all singular terms will contain a term 

Wik,u = Wtkt\ ®t k + W M t\ ®t t . 

For each of the two contributions, we can now introduce weight factors which redis- 
tribute the splitting functions to the corresponding mappings 

W lk — > A lk Wf k ] + A k iW%\ (23) 

where 

A-ik + A k i = 1 

for any fixed momenta. denotes that for the mapping of this part of the interfer- 
ence term, pi is considered to be the emitter; for W^ k \ particle p k acts as the emitter, 
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such that the roles of t and k are interchanged. We then have for the total sum of the 
two distributions 



W, 



A 



W + t 



+ A 



k(' 



wi ] t\ ® t k + t\ ® u 



We now combine this with the pure squared splitting function Wu with the colour 
factor t\ <S> tg. Invariance of the matrix element under colour rotations implies [82J: 



EH' 



t( + tn 



and the complete contribution obeying one mapping is then given by 



t e + tl ® t k 



[W u - W ik ] 



with the spin averaged interference contribution 

Pe ■ D(pj,Q) -p k 



W 



Ck 



4:7ra,2A. 



tk' 



(24) 



(25) 



Pj -PePj - Pk 

We now split the collinear and soft parts of the respective spin-averaged splitting 
functions according to 



Wu - W 



Ck 



7eik 



W u -W u ) + ( Wit ~W 



oik 



(26) 



The second part of Eqn. (|26|) can be expressed in terms of dipole partitioning functions 



4* [83J: 



W 



:ik 



w 



fk 



-P 2 

A _ A/ r tk 

4 vr a s A ek — — — — — . 

(Pj -pePj-Pk) 2 



where Py- = (pj-pe) Pk — (Pj-Pk) Pe- Several choices for A' lk have been proposed in 
all results given here have been obtained using Eqn. (7.12) therein: 



AMP} 



m+1, 



Pj -PkPe-Q 



Pj ■ Pk Pe ■ Q + Pj ■ Pe Pk ■ Q 

The partitioning weight function A' ek also obeys the relation A' ik ({p} m+ i) + 
A' ki ({p} m+ i) = 1. The general form of the interference spin-averaged splitting function 
is then given by 



AW, 



(k 



W 



■cik 



W 



tk 



Air a. 



Ipe-pkPe ■ Q 



Pe ■ Pj (Pj -PkPe ■ Q+Pe-PjPk ■ Q 
The corresponding color factor is defined by Eqn. 



(27) 



as 



Cfu = 



4 ® u + t\ ® t k 



(28) 



The only singularity in Eqn. ( 127|) arises from the factor pi ■ pj in the denominator; 
the interference term is constructed such that it vanishes for pj ■ p k — > 0. We also 
assume that the variables considered are such that they are finite for pe ■ p k — > 0, 
i.e. singularities arising in this limit should be taken care of by the definition of the jet 
function as described in Section |2] The interference term only needs to be considered if 
the emitted parton j is a gluon. If parton j is a quark or antiquark, this term vanishes. 
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3.2 Final state momentum mapping 



In this section, we will describe the momentum mapping which is used in the shower 
prescription [5TH53"] as well as the subtraction scheme. As before, hatted momenta {p n } 
are used to describe (m + l)-parton phase space and unhatted momenta {p n } m-parton 
phase space particles; emitter, emitted parton and spectator are labelled pg, pj, and pu 
respectively. The four-vectors p a ,Pb refer to initial state partons. 



For a parton splitting 

Pi -»■ Pi + Pj 

onshellness of all momenta in both m and (m + 1) phase space requires a momentum 
mapping which reduces to 

Pt = Pi + Pj 

in the singular limits; away from these kinematic regions, an additional spectator mo- 
mentum needs to be modified to guarantee pj — pj = for all particles. In our 
scheme, we use the whole remaining event as a spectator, which leads to a scaling 
behaviour ~ N 2 /2 for the number of required mappings, where N is the number of 
final state partons in the process^. 



3.2.1 Mapping in the parton shower 

For a final state splitting, we leave the momenta of the initial state partons unchanged: 

Pa = Pa, Pb = Pb- 

Let Q be the total momentum of the final state partons 

m 

Q = ^Pn = Pa+Pb- (29) 
n=l 

Here the momenta of the incoming partons remain the same, hence Q = Q = p a + Pb- 
We define 2 

at = w^-n, (30) 
2p e ■ Q 

where at > 1. The momenta of the daughter partons pe and pj are then mapped 
according to 

Pt = Pi + Pj = ^Pe H o +V Q- ( 31 ) 

The parameters A and y follow from energy momentum conservation as 

A = ^(l + y f-4a e y, y = V j^L. (32) 

y is a measure for the virtuality of the splitting, with 

, 2 



y* 



(y/a~e- y/ai-l) = 2a e - 1 - 2 y/a e {a e - 1) (33) 



3 In the Catani Seymour scheme, each additional parton in the process subsequently serves as a 
spectator, leading to an overall scaling behaviour ~ N 3 /2 for number of required mappings. 
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and A(y max ) = 0. 

The mapping prescription used in our scheme now defines the whole remaining event 
as a spectator, i.e. the momenta of all corresponding final state particles are mapped 

as 

p{J = A(K,Kr v p», n i {£,j = m + 1} (34) 
with the Lorentz transformation 

A^*, K 2 ) = ^ - ^l±Wl±Bl + ™ . (35, 



(#1 + K 2 ) 2 



.K" and _?T are given by 

K = Q-p h K = Q-P t , (36) 

and correspond to the total momentum of the final state spectators before and after 
the splitting respectively, with 

k* = A(k, Ky v k v . 

For ai — 1, this simplifies to 

k={l-y)K 

and we therefore have 

A(K,Kr(a e = 1) = {l-y)gr. (37) 
The flavours of the spectator partons remain unchanged 

fn = fn, ri£ {£,j = m + 1}, 
while the flavour of the mother parton fi obeys 

ft + fj — fa 

e.g. if the mother parton £ is a quark/antiquark, then we have (fi,fj) = (q/q,g)- 
If the mother parton £ is a gluon, then (fe, fj) can be a pair of gluons (g,g), which 
corresponds to g — > g g splitting, or any choice of quark/antiquark flavours (q, q), which 
corresponds to g — > q q splitting. 

3.2.2 Mapping in the subtraction scheme 

There is an inverse of the above mapping prescription, which maps the (m + l)-parton 
momenta to the m-parton momenta needed for the evaluation of the real-emission 
subtraction terms. We start with {p} m+ i and determine {p} m - The momentum pe of 
the mother parton follows directly from Eqn. fl3T|) 



kfo+Pi)- 1 oi +y Q- ( 38 ) 

A A A Obi 



11 



V = or> A — m and a * = op n — m » ( 39 ) 



The parameters y and read 

P 2 Q 
2P,^rpf and = 2iVQ-P- 

with = p£ +Pj. The parameter A then follows from Eqn. (132]) . 

Now we need the inverse Lorentz transformation to Eqn. (|34|) . which is used to map 
all non-emitting final state spectators. We have 

ft = h(K,ky v p v w n<£{£,j = m + l}, (40) 

where A(K, K)^ u is given by Eqn. ( 1551) . For = 1, the mapping reduces to 

Pi = (Pi -yQ),Pk= t~~ ■ (41) 
1-2/ 1 - y 

The flavour transformation is similar to the case of parton splitting. The flavour of the 
mother parton fg is given by 

ft = fi + fj , 

with the rule of adding flavours, q + g = q and q + q = g. The flavours of the spectators 
remain unchanged: 

fn = fn, n$ {£,j = m+ 1}. 
3.2.3 Phase space factorization 

In the integration of the subtraction terms over the one-parton unresolved phase space, 
we use the generic phase-space factorisation 

d{pj} m+ i g({pj}m+i) = [d{p, f} m ] d£ p g({p, f} m+1 ) , (42) 

where g({p, f} m +i) is an arbitrary function. In this work, we chose to regularise the 
infrared and collinear singularities which appear in the splitting functions using dimen- 
sional regularisation, i.e. we work in D = 4 — 2e dimensions so that the singularities 
appear as l/e 2 (soft and collinear) and l/e (soft or collinear) poles. We then have for 
the unresolved one-parton integration measure 

dt v = dy % min < y < y max ) A D ~ 3 *L<L 2 vr 5+{p 2 ) 2 vr 5+(p 2 ) 

x (2 tt) d 6™ (p t + p j -\ Pe - 1 ~ A + ?/ q] . (43) 



Here y m m = for massless partons and y max is given by Eqn. (133|) . The reduction of 
this measure for the simple case ag — 1 has been presented in [59]. In this work, we 
have used the parametrizatiorjf) 



dip 



(2 Pi -Q) 



l-e n -%+e 



16 r(|-e) 

!/max pi pi 

dyy^X 1 ' 26 dz[z(l-z)\~ e dv [v (1 - v)]~^ . 
n Jo Jo 



1 + 2 £ 



4 We thank Z. Nagy and D. Soper for useful discussions concerning the parametrization of the 
integration measure. 
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In the center-of-mass system of pi, pj, where pk defines the x — z plane, z and v 
parametrize the polar and azimuthal angle of pj respectively. 



3.3 Generalized final state subtraction terms 

In this section we present results for the subtraction terms T>e and their integrated 
counterparts Ve for final state emitters, where at ^ 1. Results for the simpler case of 
maximally two final state partons as well as initial state emitters have been presented in 
[59] . The integrated subtraction terms Ve contain integrals which depend on maximally 
two additional variables and need to be integrated numerically. In the expressions 
below, we leave out a common factor 4ira s in the expressions for the squares v\ of 
the splitting amplitudes; the (integrated) subtraction terms T>i (Vg) contain all factors. 
We will summarize the scheme in Section 13.41 We used the Mathematica package 
HypExp [86,87] in some of our calculations. 



3.3.1 Parameters 



In this paper, we use the labelling T^fjj. and Vf t f t f. for a process with the splitting 
Pe —> Pe+Pj- For final-state splittings, the subtraction terms can be expressed through 
the variables 

Pi-Pj i Pj- n £ (AA ^ 

V = and z = , 44 

Pe-Q Pe-n e 

with 

Pe = pe + Pj, n e = - Q- —P, pt ■ Q = Pe ■ Q - Pe ■ Pj , (45) 
where we additionally introduced 

l + \ + y 1 - \ + y 

7 = o , %o = 7— -7—— ■ 

2 1 + A + y 



3.3.2 Collinear subtractions 

We first consider the collinear part of the subtraction terms, which are given by the first 
term in Eqn. f[2"6"|) . These terms do not contain any soft or combined soft/ interference 
singularities, i.e. they only contain single poles ~ e" 1 and do not depend on a specific 
spectator k. 



qqg, qqg 



The squared splitting amplitude for final state qqg couplings in the case of massless 
quarks is given by 



qqg 



where 



y (pt ■ Q) 



eik 



1 



(A - l + y) 2 + 4y 
4A 

1 + Xq 



I - ~~, ~~ 2 [1 - ?y - AU ; (46) 



-1 + 



4 



Xq 



x + z(1-x ) (x q + z(1-x ))' 



(47) 
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Thus we have 



u Ana s 

QQ9 2 



C F ( 



^qqg ^cik 



C f 



(X-l + y) 2 + 4y 



y(pt-Q) r I i A 

and the integrated subtraction term is 



-Feik + 



D-2 



z[l + y + X}\, (48) 



a. 



1 1COll _ h (~1 

993 ~ 4tt F Y{\-e) \ Vl -Q 



2n fi 



2\ £ 



-+ 4/ 3 (a/) 



+ - [(9 - 7a<?)(a<> - 1) log(a^ - 1) + a^(7a^ - 16) log(a^) - 71og(y max ) 
-ae(2y max + 7) - 7y max - 4] [>, 



(49) 



with 



/•J/max 
JO 



dy 



(A-l + y) s 
4y 



+ 1 



1 + y) In 
A 



(50) 



gqq, gqq 

The gqq splitting function for massless quarks, keeping the gluon helicity for the mother 
parton, is given by 



PePj 



—Qw' — 2 



PePj 



where k± can easily be obtained from a Sudakov parametrization as 



k± = p t -j[Pn(l-z(l + x ))+yQ(2z-l)], 



(51) 



with k± ■ pi = k± ■ ni = 0. If there is no explicit helicity dependence in the Born-type 
matrix element, we have 



9qq ypeQ 

We obtain for the subtraction terms 



(1-e -2z{l-z)). 



47TO! 

(v\V gm W) = Ana s T R (u\v 2 gqq \u') : V^ gq = 



2 

gqq' 



(52) 



(53) 



Integrating this over the unresolved one-parton phase-space yields 



, x T R a s (2-k n 2 ^ ■ ' 



18 1 

+ - [(a e - 1) \n(a e - 1) - a t In a e ] 

68 9 6 



(54) 
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The total (unaveraged) splitting amplitude squared, in the helicity basis of the mother 
parton pg, is given by 



2 (Pe ■ Pj 



^— {-g^ipt ■ D j ■ p t + pj ■ D £ ■ fr) + kl ki Tr[D t ■ Dj] } , 



(55) 



with 



and 



Pt ■ Dj ■ Pe 

Pj ■ D e ■ pj 
Tr [D t ■ Dj] 

A 



2ype-Q 

x + z (1 - x ) 

2ype-Q 

1 — Z (1 — Xq 

D-2-2A + A 2 



Q 2 (pe-Pj] 



z(l -Xq) 



X + Z (1 — Xq) 



Xq 



Xq + Z (1 — Xq) _ 

Xq 

1 - z(l -x )_ 



2x 



(pi-Q)iPj-Q) ~ (xq + z(1-xo))(1-z(1-xq))' 



(56) 



and k± again given by Eqn. (15 ip: if the Born matrix element is helicity independent, 
we have 

«i» = 2( /^. )2 {P-2) ^• J D i -^ + p r ^-p i ] - klTr[D e .Dj]}, 



with 



ki 



-2yz{l-z)pi- Q. 



(57) 



Instead of using this as a subtraction term, however, we proceed in a different way, and 
define a subtraction term that only contains soft singularities from particle j [82]: We 
introduce 



(H^L,sut>'> = M v 2 ~ v lW) 



9" 



2 (p^ ■ Pj) 



[pi ■ Dj ■ p e - pj ■ D t ■ pj] 



(58) 



where v 2> 3 are defined corresponding to Eqns. (2.40)-(2.42) in [32J. This leads to 

/ l~2 I A / I 2 i 2 I a 



^— ^ { -2 <r' • Dj -p £ + kl k^Tr [D t ■ Dj] } 



2 (P/ "Pi) 

which is the subtraction term for each gluon emission. The first part is the unaveraged 
eikonal splitting function; if we combine this with the interference term, we have 



'Mv 2 



999 



rv I rv I 



2 (Pe-Pj) 



[D — 2 — A(2 — A)] . 



The collinear subtraction term reads 



2 "' n: ^7^^7 [D-2-A(2-A)]. 



l-e 



(Pe-Pj) 2 
15 



(59) 



(60) 



If there is no angular correlation in the Born- type matrix element, we can replace 
kj_ k v ± — > —k\_ in the above expressions, and equally need to multiply by 1/2(1 — e). 

Note that the above reshuffling of singular terms requires that for a final state with 
g{px)g{p2)i both combinations = (1,2), (2, 1) need to be taken into account; the 
factor | which is included in Eqn. fl55|) and all subsequent expressions guarantees a 
correct mapping of the singularity structure. 

Integrating and taking all averaging factors into account gives 
'27r// 2 \ e 1 a, 



Qj r(l-e)27r 



.4 



14 1 



(61) 



with 



hn(ai) a = X a e \ 1 - ^fa t In 1 \ _ i n 



y/at - I J \(i£ - 1 

o^=i 3 2 7 
= "gvr + (62) 

3.3.3 Soft and soft/collinear subtractions 

We now discuss the integration of the interference term, which is given by the sec- 
ond contribution in Eqn. fj26|) . This term does not depend on the specific nature of 
the splitting, i.e. it is universal; it contains all soft and soft/ colliner singularities and 
equally depends on a spectator parton k. Parton j needs to be a gluon, otherwise this 
contribution vanishes. 

We start from the definition of the interference term in Eqn. (127]) : 

1 Aw ek = Hh-h)(PfQ) 



An a 



(pt ■ Pj) [ (pj ■ Pk) (pt ■ Q) + (pt ■ Pj)(Pk ■ Q) 
and the subtraction term 

^(pi,Pi,Pk) = C lk AW ek , (63) 

where Ce k is given by Eqn. (128|) . Note that the above expression holds also for cases 
where the mother parton is a gluon, as the interference terms are diagonal in helicity 
spac^ 



3 I.e. for helicity dependent Born-type matrix elements M., where a spin correlation tensor T 1 *" is 
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We obtain for the integrated subtraction term 



with 



Vl = fi 2£ C ek I dC P (AW ( 



'ik) 



2 /x 2 7r\ £ a 



p e -Qj ir r(l-e) 



a 



ik 



i i 

K< 2? + ! 



1 + X - In (sf } + a,) 



-y + 3- 2 1n2 1n (ajf ) + 



1 

+ - 

7T 



+ In ci£ 



2 In 2 - ^ In + i In (d { Q ik) + at) + 1 | , 



(64) 



7T 

2 



1 f , 1 

— < 21n2 + 

u I ^1 + 46(1 + 6) u 2 



x In 



(1-u) 



(l + 2bu + ^1 + 46(1 + 6) u 2 )' 



+ 2 In 2 In (1 + 6) + i In 2 (1 + b) + ^Li 2 - -Li 2 

1 du f 1 dx 



b+1 



iff (««,««">) = X / £ / 
Jo " Jo 



X 



1 - X + x 0jt 



1 at 



yj [AM (1 + x 0/ -x)+ x 0/ (V^ + l)] 2 - [BW] 2 (1 + x 0/ - x) 



- 1 



+ :r - 



1 + 4. 



,-(*,fc) 



~{t,k) 
1 + ^ 



^finV*) = 71 I dx ~ ln 

JO X 



We have introduced 



at 



S e a £ 



X 



(65) 



AWfapk) = zae^r + (l-ze)~z m , B (ft) (p<,ft) = 2^^(1-^(1-^)), 



defined such that —g^jT^ = \M\ 2 , the interference subtraction term is given by 

-g^D if T^ = D if \M\ 2 . 
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and 



1 - %o,e It 

Pk -Pi Pk- Pi 

&t — \ / aj — <h — — 



St{x) = ^(l + *) 2 -l, 

2x 

where ye in all above expressions is defined b)j§ 

y e ■= S e (x)u. 

We here made the dependence on the momenta pe,p k explicit in the labelling of the 
variables A, 7, x , which are all defined according to Sections E2] and respec- 
tively. 

In terms of the Born-type kinematics, Pe can easily be recovered from Eqn. f!3ip : pk 
needs to be reconstructed in the Born-type integrations according to 

p k = A(K,K)p k , (66) 

with the Lorentz transformation defined according to Eqn. f )35|) and with 



K = Q-p e , K = Q(i-WM) - \ m . 

Finally, note that 

~(e,k) _ Pk ■ Q 



at + 



Pk - Pe 



3.4 Final expressions 

In this section, we describe how the expressions in the last subsections should be com- 
bined to provide the subtraction terms da a and their integrated counterparts j l da a- 

The complete parton level contribution is given by the sum of a^ b ° and a^ b LO , with 

<*d? = I d°ab{Pa,Pb), 
J m 

d(J ab(Pa,Pb) + / da^ b (p a ,p b )+ / da%,(p a ,Pb,fJ? F )- 

-1 J m J m 



NLO 
ab 



The NLO contribution can be split into 

[d(?ab(Pa,Pb) ~ da£ b {p a ,p b )] 



NLO 
T ab 



+ 



m+1 

d&ab(Pa,Pb) + I da^ b (p a ,p b ) + da° b (p a ,p b , 



e=0 



3 We thank Z. Nagy for providing us with this variable transformation for the interference terms. 
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where da^ b + da^ b can be written as 



dVab(Pa,Pb) + da^ b (p a ,p b ,ii 2 F 



dx / da^ b (xp a ,p b ) g) [K a (xp a ) + P(x,n 2 p 



+ dx da^ b (p a ,xp b ) ® [K b (xp b ) + P(x,n 2 F 

JO Jm 



where the insertion terms K, P only appear in the case of initial state partons, c.f. Ap- 
pendix [HI All observables, as well as infrared safety of the Born level contributions, 
need to be introduced in terms of jet functions as discussed in Section [21 c.f. Eqn. 
( II II) . For an incoming lepton, the collinear counterterm is set to zero and the PDF is 
replaced by a structure function ffft = 5(1 — r/j). 

In the following, we discuss the specific form of da^ b (p a ,p b ) which corresponds to the 
subtraction term in the real emission contribution of the process, as well as the inte- 
grated D- dimensional counterterm j 1 da^ b (p a ,p b ). In general, the subtraction term can 
be split into contributions originating from all possible emitters pt 0: 



dcr ab 



(Pa,Pb) = ^2da^ e (p a ,p b ) 



(67) 



where p( can denote an initial or final state particle. We have for each contribution 

AW 



da^(p a ,p b ) 



m+l J m+l 



»./. ./; fM'-i'i) ® \M(p)\ 



2 



(68) 



where |-M|^.y denotes the squared Born matrix element with a flavour fi of the mother 
parton; the extension for cases where there is an angular dependence of the Born-type 
matrix element is straightforward. The momenta {p m } are determined from {p m } 
through the respective mapping. N m+ i incorporates all symmetry factors of the m + l 
process and $ m +i = 2 s is the respective flux factor. For splittings where the mother 
parton is a gluon, we use the following conventions: for g — >■ qq final state splittings, we 
always choose (fi, fj) = (q, q); for g — > gg, i.e. a final state which contains g(pi)g(p2), 
we need to consider both combinations (pe,pj) = (Pi,p2), (p2,Pi)', we compensate this 
by introducing an additional factor | in the respective (integrated) subtraction terms. 
This factor has already been accounted for in all expressions in Section 13.31 



The subtraction terms can be split into collinear and interference terms: 

^r, /•,/•.(/'-•/', ) = V 7jj(P^) + &h* E ^W^), (69) 

where V li (pi,pj,pk) now denotes an interference contribution where pk acts as a spec- 
tator as discussed in Section 13.1.21 Note that there is a unique momentum mapping 

7 In the following, we omit the jet functions for notational reasons; however, full expressions should 
always be read according to Eqn. (1111) where all jet functions are included. 
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for each combination {pe,pj) which is the same for all interference terms appearing in 



The integrated counterterms are given by the integrated form of Eqn. (I67p : 

d°ab(Pa,Pb) = l^2dcr^ e (p a ,p b ). 
Ji t 

The collection of the integrated counterterms is then straightforward: for each dipole 
which has been subtracted in the real emission part, the respective integrated contribu- 
tion to /, K, P needs to be added to the virtual contribution as in Eqn. (HOOD . Finally, 
our expressions have been derived on a matrix element level: 

\M\ 2 m+1 Jv® \M\ 2 rn = V®\M\ 2 m , 

on cross section level, we additionally have to take the flux as well as combinatoric 
factors into account)! 

dai +1;a ,(p a ,p b ) = fv ® \M\l = ^g±V ® \M\l, 

dv m+1 -, ab = N m+1 [ iv ® \M\ 2 m = %±i* s V® / da m , 

where the factors N m , N m+ i account for possible symmetry factors of the specific pro- 
cess, and where here x s = s/s is the ration of the partonic center-of-mass energies 
before and after the splitting; x s = 1 for final state emitters. We then obtain the 
relation 

£V= + K + P) (70) 

between the integrated splitting functions V given in the next sections and the insertion 
operators J, K, P. 




4 Example: e + e — > 3 jets 

In this section we consider the simplest non-trivial process with more than two par- 
tons in the final state: three-jet production in e + e~ annihilation. The next-to-leading 
order contributions to this process are well-known [65 tlg5H5U] . We compare the results 
obtained from the implementation of our scheme and from a private implementation 
of the Catani Seymour scheme as well as [9oJE We find complete agreement for the 
differential C-parameter [63], with integration errors on the percent level. 

The leading order process we consider is given by 

e+e" — > q(Pi)q(P2)g(P3)- (71) 

8 Correct counting of symmetry factors needs to be done explicitely in this expression; if all splitting 
multiplicities and symmetry factors are taken into account, we obtain a generic combinatoric factor 
\ for ggg splittings, c.f. Section 7.2 in [55] , 

9 We thank M. Seymour for help with the original code available from [91 . 
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At next-to-leading order, two different real-radiation subprocesses contribute: 

(A) e + e~ -> 7*(Q) -> g(ft) g(p 2 ) g(p 3 ) g(p±), 

(B) e + e~ 7*(Q) -»■ g(ft) g(p 2 ) g(p 3 ) q(Pa)- 



(72) 



For a complete next-to-leading order calculation, the virtual corrections need to be 
added to the leading order contribution. In the following, we use the notation 

2 ft Q / Q , t/ij Q j Sijk &ij ~\~ Sjfc -|- Sj] il 

with Sij = (pi +pj) 2 . Energy- momentum conservation leads to 

J^Vij = 1, J2 Xi = 2. 

i,j>i i 



We equally follow the notation for matrix elements in Section 13.11 



( M I M ) 

(m\ {Pi}) 



(M 3 ({ Pt })\ Ms({ Pi })) = \M a (fa})f 



(Mt{{pi})\ M,({pi})) = \M,({p z })\ 2 
The total next-to-leading order contribution for process (171 j) is then given by: 



a 



NLO 

J 



dPS d 



+ J dPs A 



flavours 



\M* ({Pi})\ 2 Ff ({ft}) - (Q, r, s\ V £j ({ft}) \q, r, s) F ( f\q, r, s 

e,j 

\M 4 ({pi})\ 2 Ff ] ({p t }) -J2(q,r,s\ V ej ({ft}) \q, r, s) F®{q, r, s 
+ / dPsJ \M V ({ Pl })\ 2 + ^ \ ^(1,2,31^- ({ft}) |1, 2, 3) 

+ E i E^/dep<i,2,3|^(^ i })|i,2 ) 3> Uf ({ft}) 

flavours L ' J p J 



(73) 



where £,j sum over all possible pairs in the real emission phase space, ft and ft denote 
momenta belonging to Born- and real emission kinematics, and where q, r, s G {1, 2, 3} 
denote the mapped momenta in the real emission subtraction terms for the Born-type 
matrix elements. The factors ~, \ in the integrated subtraction terms correspond to the 
process-dependent symmetry factors in the real emission contributions. In the symbolic 
notation above, J dPS n contains all symmetry and flux factors for the respective phase 
spao 



4.1 Tree level result 

In the following, we follow the procedure of 
according to 

1 

— do j 



i.e. we normalize our observables 



da^ LO 



10 



Note that / dPS 3 , J dPS A are defined slightly differently in [551150] . 
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(74) 



where <r denotes the total cross section for the process 

e + e~ — > q q 
given by [63J 

Co = — — L>Aq f 

o S 

for a quark with flavour / and charge qj. We normalize the matrix elements according 
to 



a 



(n) 



J dT n \M n \ 2 F 



2 jp(n) 
J > 



where dT r 



mass energy and with 



0fz 5 (pf - mf) 5 (4) (J2 in Pm ~ J2i Pi), s denotes the center-of- 



S 



1/2! for 7* — ► qqgg 

1/2! x 1/2! for 7*-»gggg 



being the symmetry factor of the respective process. We obtain the well-known rela- 
tions between the matrix elements of the processes (1711) and (|71|) [S3] : 



C h 



ret •-' I rf*£- 
dj 1 \ JL O 



1 



[l-x 1 )(l-x 2 ) 4 



where \Ai 2 1 2 has been averaged over the emission angles, as well as 

o"3,j = tt-Cf / dxidx 2 Ih [0(1 - Xi)Q(xi)] Q(x 1 +x 2 -l) , _ 1 v 2 _ — -a F}\x 1 ,x 2 ) 

Z 7T J yL XijyL X 2 ) 

for jet observables. The gluon-helicity dependent squared matrix element for the lead- 
ing order process (TTTj) is given by [90J 



where 



r^(pi,p 2 ,p 3 ) 



--T—r, \M 3 ( Pl , P2 , P3 )\ 2 T^, 



rpfiv = +2 PlPl + 2 _ 2 1 "Si PlPl _ 2 1 ~ X 2 



g 2 



Q 2 



1 — X\ — x 2 + x\ 
\ — x 2 



1-X2 g 2 

' A 1 f A* ^ 

, g 2 g 2 



l-xi g 2 

1 — x 2 — Xi + x\ 

\ — X\ 



' A 1 v V- v 

P2P3 , 

. g 2 g 2 



(75) 



(76) 



Tut, obeys g^T 



111/ 
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The virtual one-loop matrix element 

The virtual contribution to the process (I7ip in the MS renormalization scheme is [B^|Ido] 



I a I \i2 i a I \i2 a s ( 4VTyU \ 1 

\Mv{Pl,P2,P3)\ = \M 3 {Pl,P2,P3)\ 77- 



2tt V Q 2 J r(l-e) 



X 



l -i [(2C F - C A ) y£ + C A (y^ + ^)] - i (^3 C F + ^C A - |n, T R ^j 



7r 2 „ s I a. 



+ - (2 C F + CU) - 8 C F \ + [F(y 12 , y l3 , y 23 ) + 0(e)] , (77) 

L J Z7T 

where -F (2/12, 2/13, 2/23) is defined in Eqn. (2.21) of [63J. 



4.2 Subtraction terms 

In general, the subtraction terms are given by Eqn. f lo^]) as 

where the first part is the collinear subtraction term, and the second the sum over all 
interference terms. For the processes considered here, the interference terms only con- 
tribute to (A), while (B) only contains collinear divergences. In the following sections 
we will construct the subtraction terms for subprocesses (A) and (B), respectively. 



4.2.1 Momentum mapping 

For all subtraction terms, the m-parton phase space momenta are mapped as described 
in section l3.2.2( i.e. , for a splitting 

Pe ->■ Pe + Pj, 

we have 

Pe = -rPe 7T\ Q> 

A 2 A dp 

Pk = A(K,K)p k , [k (£,j)] 
where A is defined according to Eqn. ([35]) . and where 

K = Q-Pt, K = Q- Pt . 



4.2.2 The subprocess: 7* ->• q(p 1 ) q(p 2 ) g(p 3 ) g(p 4 ) 

For this process, ggg, qqg collinear as well as interference terms need to be taken into 
account; these are given by Eqns. (pi]) . (SSJ), and (153]) respectively. We have for the 
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collinear subtraction terms 

4-7T a 



[v\V 



Coll I , ,/\ T-I/l/ 



999'M 



|^ 3 (Pl,P2,P3)| X 
1 — ^2 



^ ( 2(*jl • Pl )(A; ± • p 2 ) - • Pl ) 2 - —^(k ± ■ p 2 ) 



+ k l ( 1 + y + y - *i " ^2 



with fc^ and A given by Eqns. (I5ip and ( |56|) . In total, we need to consider the following 
combinations: 

<T-)C0ll <T-\C0ll -T-JCOII -T-JCOH <J-\Coll <T-\Coll 

U qqg;\Zi ^ qqgMi U qqg;23> U qqg:2^ U ggg;3Ai ^ggg;43> 

note especially that both combinations (3,4), (4,3) need to be taken into account in 
the contributions from ggg splittings. 



The interference terms are given by 

Vf J;k =AW £j (k) = An a a a 



2 P i -PkPe ■ Q 



£k 



Pi ■ Pj [Pj -PkPi- Q +Pe- Pj Pk ■ Q 



(7c 



where we have the following contributions: 

7}if W -7-)if r>if T)if r>if -T)if r>if r>if T)if r>if T)if 

■^1,3;2' ^1,3;4> ^1,4;2> ^1,4;3' ^2,3;1> ^2,3;4' ^2,4;D ^2,4;3' ^3,4;2' ^4,3;1' ^4,3;2- 

We want to emphasize again that there is only one mapping required for each (£,j) 
pairing, i.e. we only have 5 independent mappings for the 12 subtraction terms listed 
above 11 !. The subtracted cross section is then given by 



a R ~ A = J [da R - da A ] 

= J dPS A | |A^ 4 (Pl,P2,P3,P4)| 2 F 4 (p 1 ,p 2 ,P3,P4) 

-{ 1, 2, 3 I Dis + V U + £>23 + £>24 + £>34 + £>43 1 1,2,3 )F 3 ( Pl ,p 2 , p 3 ) } 

(79) 

with the spin-averaged matrix element 

\ \M 4 (pi,P2iP3iP4)\ 2 = CTO dp A (4 7T) 5 X 

{(A + 5 + C) + (1h2) + (3h4) + (1h2,3h4)} (80) 

The quantities A, B, C are given in Appendix B of [63] . 

11 The Catani Seymour prescription (65) requires 10 different mappings for this contribution. 
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4.2.3 The subprocess: 7* — > q(pi) q(p 2 ) q(Ps) q(p4) 



This process does not contain any soft/ interference singularities, and all subtraction 
terms are given by Eqn. (I53p . We need to keep track of the mother partons helicity in 
the subtraction term and obtain 



[v\V mq ,^\v')T vv ' = 4<rca s Z 



R 



|-M3(Pl,P2,P 3 )| 2 X 



Pi ■ Pj 



{k± ■ P2) 



+ k l ( 1 + y + y - x i ~ x 2 



with k± and k\ given by Eqns. ( 151]) and (|57|) respectively. 
We have to consider the following combinations 



gqq;l2i T^gqq-M, ^ 'gqq;32, ^gqq;3A 



(81) 



and get 



R-A 



= j [da R - da A ] 

dPS 4 j \M 4 {pi,p2,P^P4)\ 2 Fi(Pl,P2,P3,n) 

-(l,2,3|Pi3 + V M + V 32 + V 3A \ 1, 2, 3 )F 3 (p 1 ,p 2 ,p 3 ) } 



with the spin averaged matrix element given by 

1 l^4(pi,P2,P3,p 4 )| 2 = ^0 (^) C F A (4vr) 5 x 

{(£> + i?) + (1 fy 2) + (3 H 4) + (1 « 2, 3 f> 4)} 

The quantities D and E are given in Appendix B of [63J. 

4.3 Integrated subtraction terms 
4.3.1 Collinear integrals 

The collinear integrals involve the gqq, qqg/qqg and ggg splittings: 

^coii = i co \\{gqq) + i co n{qqg) + honiggg)- 



(82) 



54) 
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After summing over all contributions, we obtain from Eqns. (154"]) . (l4"9]) and (16T1) 

(l,2,3|I coll (<?^)|l,2,3) B = 



\M s ( Pl ,p 2 ,p 3 )\ 2 ^T R 



2\ e 



1 



ra-el 



71/ x 



2 16 2 

•- — + - [(a 3 - 1) ln(a 3 - 1) - (a 3 + 1) In a 3 ] 

o s y o 



(1,2, 3 1 10^^)11,2, 3) A = 



2\ £ 



ra-e 



x 



^ ^ < - - + 4 1 3 (d£) - In + ^ [(9 - 7a e )(a e - 1) ln(a £ - 1) + a e (7a t - 16) log(o^) 
-7log(y max (at)) - a e (2y nmx (a e ) + 7) - 7y max (a £ ) - 4] 



;i,2,3|/ coll (^(?)|l,2,3) A ) = 

14 1 1 

-7 n + 7 t( a 3 - 1) ln («3 - 1) - (a 3 + 1) In 03] + / fln (a 3 ) ^ , ( 

6e 9 6 J 

where all symmetry factors have already been taken into account. Here, and 7 3 are 
given by Eqns. (16"2"j) and (150)) respectively and need to be evaluated numerically. 



55) 



4.3.2 Soft integrals 

For a specific emitter/ spectator pair (pe,Pk), we obtain from Eqn. (IMj) 
4 7r a, 



'soft,flb 



a s ( 4:77/1 



2\ e 



ra-e 



2e 



1 1 
+ - 



1 + - In 

2 



7T 



+ 3 



+ 



7T 



'fin 



+ 41 ) (^ ) ,%)+4 e n ) (^) 



+ 7 ln ' 
4 



~(ik) 

- In 2 In I Q ° I + 2 In 



(86) 



where 



(6)/ 



In 2 ln(l + 6) + - ln 2 (l + b) 



The additional factor | in the integrated subtraction terms stems from Eqn. (1701) 
and accounts for the different symmetry factors and combinatorics of process (A). Soft 
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interference terms only appear for this process, where the following emitter/ spectator 
pairs need to be taken into account 

2 [(1, 2) + (1, 3) + (2, 1) + (2, 3) + (3, 1) + (3, 2)] ; (87) 

the factor 2 arises as e.g. (pe,pj,pk) = [(P1P3, P2), (pip4, P2)} are mapped to the same 
Born-type kinematics (pe,Pk) = (Pi,P2), an d similar relations hold for the other con- 
tributions. 



4.3.3 Finite parts 

Combining the one-loop matrix element Eqn. ( I77p with the integrated subtraction 
terms, all poles in e cancel. The leftover finite parts are 



9 Qt 

< 1,2,3 |/ coU | 1,2,3 > finitc = \M 3 {pi,p 2 ,p 3 )\ ^ x 



(4n f T R + C A ) 



~t + l [ ^ 3 ~~ X) ln(fl3 ~~ X) ~ (fl3 + X) ln as] 



+ C A I &n (a 3 ) 



£=1,2 . 



4 1 3 (ae) - ln a £ + - [(9 - 7a e )(a e - 1) ln(a £ - 1) + a e (7a e - 16) log(a^) 



-71og(y max (a^)) - at{2y max (ai) + 7) - 7y max (a^) - 4] 



and 



< 1,2, 3 I /soft 1 1,2,3 > finite = \M 3 {p 1 ,p2,p 3 )\ 2 7T- x 

Z IT 



(2C F + C A ) 
-2(C A -2C F ) 
-I 



7T 

- T + 6 



+ 2 -£ c? 

£=1,2,3 



-I^(a e ) + 2\n a e 

7T 



1 =f. 



~(1,2) 
•- j(b) [ % 

fin 1 

7T \ di 



In 2 ln 



a 1 



+ 1 



+ 2CU ]T 



=1,2 



(6) / «0 



_ fin . 
7T \ 



In 2 In I — 



(£,fc) 



(a[f' fc) + ciij a, 



39) 



where the sum in the last line goes over all possible combinations as given in Eqn. (J? 
(the factor 2 is already accounted for) , and where we made use of several symmetries. 12 



12 One useful relation is e.g. 



~(m,n) ~(n,m) 
% _ % 
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Further simplifications for the interference term finally render 



< 1,2,3 | /soft | 1,2,3 >finite= \M3{pi,p 2 ,P3)\ 2 ^ X 



(2C F + C A ) 



+ 2 Y. C l 



e=i,2,3 



In ai In a 2 — In a 3 In 



ai a 2 



a 3 



-l£>(a e ) + 2 In a e 

71 



+ In 2 ?/i3 + In 2 ?/23 - In 2 yn 



+ Ci 



2 

+ - 

7T 



a 2 



ln 2 (^) + ln 2 y 12 



+ 2 In 2 



2 Ci? In (yu ai a 2 ) - CU In 



1/12 



yi3y23«3 



fin 



=1,2 



1 \ 



(90) 



4.3.4 Final expressions 

The finite parts of one-loop matrix element are given by 

2 2 
\Mvipi,P2,P3)\ &nite = \M3{p 1 ,p 2 ,P3)\ —X 



| X - [ (C A - 2 C F ) In 2 j/ 12 - C A (In 2 y 13 + In 2 y 23 ) ] + y (2 C F + CU) - 8 C F J 



+ 



2tt 



F(y 12 ,y 13, y2s) + 0(e) 



(91) 



If we combine the one-loop matrix element with the integrated subtraction terms, poles 
in e exactly cancel, leading to finite results: 



r V+A 



d(T V + J^da' 



1 



= / dPS 3 { (l,2,3|/ soft |l,2,3) + (l,2,3|/ coll |l,2,3) + |7Wy(p 1 ,p 2 ,p 3 )| 2 }F 3 (pi,P2,P3) 



= / dPS 3 ^ (1,2,3 |7 soft | 1,2, 3) fin + (1,2,3 |/ coll | 1,2, 3) fin + \M v (pi,P2,P3)\L }F 3 ( Pl ,p 2 , 



(92) 
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where 



(1, 2, 3 1/^ 1 1, 2, 3) fin + (1, 2, 3 |/ coll | 1, 2, 3) fin + \M v ( Pl ,p 2 ,p 3 )\ 



fin 



\M 3 {P1,P2,P3)\ 



2 

2tt 



X 



7r 2 ^0 16 7 

(2CV + CU)— + 2C F + —C A - -n f T R +-C F \n( ai a 2 ) 



23 2 

— C A - -n F T R )\n a 3 + (An f T R + C Aj 



- [(a 3 - 1) ln(a 3 - 1) - a 3 In a 3 ] 



+ c A h n (a3)+ -if Yl 



= 1,2 



1 



4 I 3 (a^) + - [(9 - 7a £ )(a£ - 1) ln(a £ - 1) + a^(7a^ - 16) log(a^) 



-7\og(y max (a e )) - a^(2j/ max (a^) + 7) - 7y max (a^)] 



In cii In a2 — In a 3 In 



«3 



C F In 2 



+ 2 In 2 



2/12 



2 C F In (2/12 ai a 2 ) - CU In 

2/131/23 a 3 



c, E 4? ( 

4=1,2 ^ 



2 

+ - 

7T 

+ 7r- ^(2/12,2/13,2/23) • 

2 7T 



1 \ 



2te^a 3 



(CU-2CV)/, 



(6) 
fin 



£=1,2,3 
1 



2/120102 



(93) 



4.4 Results 

We compared the implementation of our scheme with the results from [65|l90] as well as 
our own implementation of the Catani Seymour scheme. The subtraction terms for the 
latter are well-known and will not be repeated here. In order to fulfill the jet-function 
requirements in Eqn. ffTUl) . we chose the variable 



1 - 



E 



i,j = l,i<j 



(2 Pl -Q) (2p r Q) 



(94) 



which is infrared finite as required_£|. We numerically compared all different color 
contributions N C C F , N c C F nfT R , N^C F separately, as well as the combined contri- 
butions. We set rif = 5 in our calculations. Figures [3] to [5] show that we agree with 
results obtained using the Catani Seymour subtraction scheme on percent-level, and 
are consistent with 0, within the errorbars, and thereby successfully validated our real 
emission subtraction terms as well as all interated counterterms proposed in this paper. 

Results for integration as well as differential distributions have been obtained using 
routines from the Cuba library [93J. 



A closer inspection of this variable shows that it contains singular regions which however are 
integrable [92] : we thank B. Webber for pointing this out. 
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C - distribution: ( C / a > da / dC 



C - distribution: ( C I a a ) da / dC. single contributions 



1000 

> 900 

^ 800 
£ 700 

a 600 

| 500 
O 400 

"5 200 
D 

O iOO 




B -400 
-600 





NS dipoles 1 — 1 — 1 




CS dipoles ' — ! - 




CS dipoles, event2_3.f 











NS dipoles i- 
CS dipoles L - 
CS dipoles, event2_3. f ' -- 



NS dipoles ' — 
CS dipoles L - 
CS dipoles, eveni2_3.f '- 



0.1 0.2 03 0.4 0.5 0.6 0.7 0.8 0.9 
C 



0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
C 

C - distribution: ( C 7a ) da 7 dC 



0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
C 

C - distribution: (C7a )da7dC 



_ 005 r 

§ 0.004 - 

v 0.00J 

^ 0.002 

a o.ooi 

a 



CJ 

5 -o.ooi 

b 

2. -0.002 

£ -0.003 
-0.004 



0.00J 
? 0.002 

5 o.oo; 
1 

-~ -0.001 

CJ 
13 

O -0.002 - 

-5 -0.003 - 
CJ 

-0.004 - 



0.7 0.2 0.3 0.4 0.J 0.6 0.7 0.8 0.9 1 

c 

C - distribution: I C / a ) da / dC, single contributions 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
C 

C - distribution: I C / a ) da / dC, single contributions 



N 0.004 - 




Figure 3: From top to bottom: Comparison of the contribution proportional to 
Nc CpTifTji, NqCf-i Cp for our scheme (NS) and Catani Seymour dipoles (CS). 

Shown is the differential distribution ^ dcJ dc° m un ^ s °f (fj) > with C defined by 
Eqn.f l94p . Left: Results for the implementation of our scheme and Catani Seymour 
subtraction terms from our private code as well as event2_3.f |91j. Right: relative dif- 
ference between our implementation of the CS and NS subtraction terms. The results 
agree on the percent level and are consistent with zero within the integration errors. 
Large errors arize in regions where the absolute values of the differential distribution 
become small. 
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Figure 4: Total result for differential distribution y- da ^° using both NS (red) and 
CS (green) dipoles. The standard literature result obtained using the CS scheme is 
completely reproduced with the NS dipoles. 
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Figure 5: Differences Aqs-ns for real emission (red, upper) and virtual (green, lower) 
virtual contributions, showing that especially for low C values the contributions in the 
two schemes significantly differ. Adding up A real + A virt gives as expected. 
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5 Conclusion and Outlook 



In this work, we have proposed an alternative next-to-leading order subtraction scheme, 
where the subtraction terms have been derived in the context of an improved parton 
shower. In this scheme, we employ a momentum mapping which reduces the num- 
ber of reevaluations of the underlying Born matrix element with respect to standard 
schemes [651 [70] 0- Furthermore, the use of subtraction terms based on the shower 
splitting functions promises to facilitate the matching of parton-level NLO-corrections 
with the improved parton shower. While final state subtraction terms for processes 
with maximally two partons in the final state as well as initial state subtractions were 
already presented in previous work [58H60] . we here extend the scheme to processes 
with an arbitrary number of final state particles and provide formulae for the corre- 
sponding final state subtraction terms and their integrated counterparts. We validate 
our expressions by reproducing the literature results for the differential distribution of 
the C-parameter at NLO for three jet production at lepton colliders, where we find 
numerical agreement between results from the implementation of our scheme and two 
independent implementations of the Catani Seymour scheme on the (sub)percent level. 
Combining the results in this work with the discussion in 15814611] provides all formulae 
needed for a generic application of our scheme for massless emitters, and therefore 
concludes the discussion of the subtraction scheme in the massless case. 

As argued in the introduction, subtraction schemes can generically differ in the non- 
singular structure of the dipole subtraction terms as well as the mapping between real 
emission and Born-type kinematics, which guarantees onshellness and energy momen- 
tum conservation in both phase spaces. The scheme proposed here uses the whole 
remaining event as a spectator for the mapping, thereby leading to a scaling behaviour 
of Born reevaluations ~ N 2 /2, where N is the number of final state partons. However, 
this simplified mapping equally induces integrated subtraction terms with finite parts 
which exhibit a non-trivial dependence on the integration parameters of the unresolved 
one-parton phase space. In this work, we chose to evaluate these finite terms numeri- 
cally, which leads to an increase of integration variables by two in the numerical imple- 
mentation of the scheme. However, recently it was proposed [lj to approximate similar 
finite terms by polynomial functions in the context of an next-to-next-to-leading order 
subtraction scheme [9H - I97] . We therefore plan to make the finite remainders which 
appear in the integrated subtraction terms available either in form of approximating 
functions or a librarized grid interpolating between different input parameters for a,£. 
Further plans for future work include the extension to the massive scheme as well as 
the matching with the improved parton showeiP^I. 



14 We want to note that the FKS subtraction scheme [53] exhibits a scaling behaviour similar to 
our scheme. However, the two prescriptions differ in the treatment of phase space setup; in addition, 
no parton shower proposal exists using FKS splitting functions. We thank R. Frederix for helpful 
discussions regarding this point. 

15 We note that work is underway [55] to implement the scheme presented here into the Helac NLO 
Event Generator framework |15j . We thank the authors for bringing this to our attention. 
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£ ft fi fj Vi x 7 === colour 



1 



... Au ^(P* + fijtyeUipe, s t ) 
F q q g £n(Pj,Sj;Q) ' , 21 t a 

[pi + PjY 



Table 1: Splitting amplitudes v#({p, f}j, Sj, se, se) involving a qqg splitting. We have re- 
moved a common factor ^Aira s . The label t denotes final-state indices F = {1, . . . , m}. 
The light-like vector ne is defined in Eqn. (I97p . Taken from 



A Splitting amplitudes 

In Table [U we list the splitting amplitudes for final state qqg splittings as given in 
For triple gluon splittings, we have for the final state 

i>t({p, f} m +l,Sj,S e ,Si) 



= 47Tas e a (p j} sj; Q)*ep(p e , s t ; Q)*e v (p e , s e ; Q) (95) 
^Pj 'Pi 

x v^ipjipt, -pj - p e ) D JV (pi + pj; m) . 

We use standard notation where U(p, s), U(p, s), V(p, s), V(p, s) denote spinors of the 
fermions with a four-momentum p and spin s, and e a (p, s; Q) are the gluon polarisation 
vectors. The ggg vertex has the form 

V a ^( P a,Pb,Pc) = 9 a \Pa ~ Pbf + 9 M {Pb ~ PcT + 9 1a (Pc ~ Paf ■ (96) 

The transverse projection tensor D^ u (pi —pj] rig) is defined according to Eqn. (122]) . The 
light-like vector ri£ is given by 

n z = Q-7\ Q2 f77^2 Pi >^{l,...,m} . (97) 

Q-Pe+ ViQ-Pe) 2 

B Incoming hadrons 

In case of processes with two initial-state hadrons A and B carrying momenta pa and 
PBi respectively, the calculation of the QCD cross sections must be convoluted with 
parton distribution functions fi/i(rji, jJ? F ) which depends on the factorization scale y,p'. 

<?{PA, Pb) = Y] / dr] a f a/A (Va, /4) / d Vb fb/B(Vb, /4) [VabiPa, Pb) + ^°{Pa, Pb, /4 
a,b J ° J ° 

(91 
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where p a = r] a pA and p b = "i] b p B are parton momenta, while r] a and 7] b are the 
momentum fractions of the partons. In this case, additional collinear counterterms need 
to be added to the integrated subtraction terms^], and parton level NLO contribution 
becomes 



We then have 



d(?ab(Pa,Pb) + I do v ab (p a ,p b ) + / da^ b (p a ,p b ,^ F ). (99) 



m+l 



Vab L °(Pa> Pb, /4) = / [ dcr ab(Pa, Pb) ~ dcr£ b ( Pa , Pb. 
J m+l 

d(?ab(Pa,Pb) + J^d(Tab(Pa,Pb) + ^S(Pa, Pb, /4) 



(100) 



with 



d&ab(Pa, Pb) + dcTa b (p a ,Pb, ^4) 



+ 



d<Ja b {PaiPb) ® 1(e) + dx da^ b (xp a ,p b ) <g> [K a {xp a ) + P( 

^ , f^F 

JO Jm 

dx / da^ b (p a ,xp b ) g) [K b (xp 6 ) + P(x,/i|)] . 



(101) 



This equation defines the insertion operators 1(e), K(x), P(x;jj? F ) on an integrated 
cross section level. Eqn. fllOip can be divided into two parts: the first part is the 
universal insertion operator 1(e), which contains the complete singularity structure of 
the virtual contribution and has LO kinematics. The second part consists of the finite 
pieces that are left over after absorbing the initial-state collinear singularities into a 
redefinition of the parton distribution functions at NLO. It involves an additional one 
dimensional integration over the momentum fraction x of an incoming parton with the 
LO cross sections and the x-dependent structure functions. 



In the MS scheme, the collinear counterterms are given by 

/ da^ b (p a ,p b ,fi 2 F ) = ^ — rV / dx I da^ b (xp a ,p b )- 
Jm 2tt r(l -£) ^ J J m e 

7 i~7) ^2j dx J daZ(p a ,xp b )^ ^ 



47T/Z 

/4 



2\ e 



P ac (x) 



2tt n 



4?r// 



2\ £ 



P bc (x). 



(102) 



Here the P ab (x) are the Altarelli-Parisi kernels in four dimensions |78] , which are evo- 
lution kernels of the DGLAP equation [75 |ll00lTl02] . and describe the behaviour of 
parton splittings by giving the probability of finding a parton of type b with momen- 
tum fraction x in a parton of type a in the collinear limit: 



16 



a(p) 

a. e.g. eng. 



b(xp + k ± + 0(k\)) + c ((1- z)p-k ± + 0(k 2 ± )) 



(103) 
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At leading order, the splitting functions are given by 
pii( x ) = CV 



1 + x 2 3.,, 



P 9g (x) 
P q9 {x) 

P"(x) 



T R [x 2 + (1 - x) 2 ] , T R 
l + (l-x) 2 ' 



1 

2' 



C F 
2C A 



x 
x 



(1-x). 



+ h x (1 — x) 



+ 6{l-x] 



llCA-AnfT, 



R 



6 



(104) 



where n/ is the number of quark flavours in the theory. The + distribution is defined 
in the standard way 

/ f(x)g+(x)dx = / g(x) (f(x) - /(l)) dx = / g(x) f(x) dx - f(l) / g(x) dx 
Jo Jo Jo Jo 

(105) 

for the convolution with a test function f(x). 



C Four-particle phase space 



In this section, we derive the parametrization which was used for the real emission 
phase space in Section |H We use the standard notation for an n-parton phase space 
in four dimensions: 



dY r 



n 



(2vr)< 



5 (4) 



Ph 



Pi 



We build our parametrization from a successive chain with 

Pin -> Pl2 + P34, Pij -> Pi+ Pj, 

where in the first step the onshell-condition for p^ needs to be replaced by a distribution 
Of Sij = (pi +Pj) 2 . 



Generic massive two parton phase space, center-of-mass system 

For a generic massive two parton phase space, we use the following parametrization in 
the center-of-mass system 



dT, 



1 y/\(s,ml,m$) 



32 7T 2 



dVLi (-y/i - (mi + m 2 )) 



where 



s + m 2 — m 2 , JA(s,m?,mo) 

Pi = ;r~F ,P2=Pin-pi- 



2^s 7 1 J " 2^ 
The 6 function arizes from the conditions 0(| p i|) 0(|"~z7" 2 |)- A is defined as 

x(x 1 ,x 2 ,x 3 ) = x ■ - 2 . 



3/ 2 • 
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Generic massless two parton phase space, non-center-of-mass system 

If we consider two partons in a non center-of-mass system, with 





Ph 







V l^inl / 

i.e. where the three vector of p in determines the z-axis, we obtain for a two parton 
massless phase space 

I 

16 7T 2 \p in \ 

where we have 

If m 2 \ 
cosO, = Ie*--*) (106) 

I P m| \ ^ Pi J 

and 



mi 



2(£ in +|^ in |) " Pl " 2(£ in -|^m|) 
from the requirement that | cos#i| < 1. In this derivation, 

Pi, z = Pi cos O-i. 

If the z-component of p m goes into the negative ^-direction, cos 9 — > — cos 9 in Eqn. 
f ll06p . and all other above relations still hold. 



Generic four parton phase space with massless final states 

We use the generic expression 

/ \ ^ \ dm\ dm\ 

dT n —> 2_^Pn) = dTx^Y+z — ^— dT X ->j2Po « r y_>£x, 

where J2p + J2p'o = J2Pn is the sum over all n outgoing particles. Using the 
expressions above as well as 

Sij „Pi ' Q 

s s 

we obtain for the four-parton phase space in the center-of-mass system of p m = Q: 

dy 12 dy^ dxi dx 3 d<f>s, 



<fT 4 

with the four-vectors 
P12 = 



(4^)6^(1,^12^34) 



P34 



(E X \ 



\ Px J 

( 1 \ 



/ Ey \ 




V -Px ) 



Pi = xi 



sin 9i 


y cos^i i 

P2 = Pl2 ~ Pi, Pa = £>34 - P3 



( 



P3 = X 3 



\ 



sin #3 cos 03 
sin 6*3 sin 3 
V cos 6» 3 J 
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and where 



P y/s n . x v/sA(l, 1/12,1/34) 
£x = -y- (1 + yn ~ 2/34) , Px = 2 ' 

= y/s- E X = — (1 + 2/34 - Vu) , 

ft 1 m x\ 1 A , 2y 12 

COS 0! = £ X p = - == ( 1 + J/12 — 2/34 - 



V X! VsJ y/\(l,yi2,yu) V X! 

COS 6*3 = £y p = = 1 + I/34 - 1/12 



px V ^3 vv v/A(i, 1/12, 1/34) V ^3 

The integration boundaries are given by 
2/34 < (1 - v 7 ^) 2 ' 



min /max 

m x 2 1/12 



•^3 



v/s (E X ± Px) 1 - I/12 + I/34 ± a/A(1, I/12, I/34)' 

min/max m Y 2 1/34 

V^^yipx) I-I/34 + 2/12 ± v /A ( 1 ? 2/i2, 2/34)' 
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